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A NEW DATA-PROCESSING TECHNIQUE FOR THE ELIMINATION 
OF GHOST ARRIVALS ON REFLECTION SEISMOGRAMSf 



WILLIAM A. SCHNEIDER,* KENNETH L. LARNER,: J. P. BURG,* and 
MILO M. BACKUS* 

A new data-processing technique is presented for the separation of initially up-travciing (ghost) energy from 
initially down-ira vcling (primary) energy on reflection sctsmogr&ms. The method combines records from two or 
more shot depths after prefiltering each record with a different filter. The filters are designed on a least mean squarc- 
crror criterion to extract primary reflections in the presence of ghost reflections and random noise. Filter design is 
dependent only on the difference in uphote time between shots, and is inde(>endent of the details of near-surface 
layering. The method achieves wide band separation of primary and ghost energy, which results in 10-15 db greater 
attenuation of ghost reflections than can be achieved with conventional two- or three-shot stacking (no prefiltering) 
for ghost elimination. 

The technique is illustrated in terms of both synthetic and field examples. The deghosted field data are used to 
study the near-surface reflection response by computing the optimum linear filter to transform the deghosted trace 
back into the original ghosted trace. The impulse respKinse of this filler embodies the effects of ihe near-surface un 
the reflection scismogram, i.e. the cause of tnc ghosting. Analysis of these fillers reveals that the ghosting mecha- 
nism in the field test area consists of both a surface- and base-of-wcathcring layer reflector. 



INTRODUCTION 

In exploration seismology when an explosive 
source is detonated below the free surface, the 
initial up-traveling energy is either absorbed or 
returned to the subsurface by the near-surface 
reflection complex. If the reflected energy forms 
a coherent wavefront, then we may expect to 
record both primary and ghost arrivals from 
every reflector at depth, separated in time by ap- 
proximately twice the travekime from the shot 
to the near-surface reflection complex. This 
"double image" view of the subsurface which 
results from ghosting is clearly undesirable for 
the following reasons: 

(1) Doubling the number of reflection events 
on the record compounds the subsurface interpre- 
tation. 

(2) Ghosts and primaries may interfere de- 
structively, thereby masking subsurface informa- 
tion. 

(3) The ghost-reflection interference may vary 
considerably within a prospect, causing poor char- 
acter correlation. 



The latter is perhaps the most undesirable fea- 
ture of ghosting. Since the ghost reflections are 
sensitive to near-surface conditions in the vicinity 
of the shotpoint, such as thickness of the weather- 
ing layer, roughness of surface and base-of- 
weathering reflector, depth of the shot, position 
of water table, to mention a few, we might well 
expect ghosting to vary significantly within short 
lateral distances. This variation in the ghosting 
can destroy meaningful fine structure in the pri- 
mary reflections which may convey stratigraphic 
trap information. For these as well as other rea- 
sons, ghost arrivals are classified as seismic noise 
along with surface waves, multiple reflections, 
scattered energy, wind noise, etc. 

The elimination of ghost energy from seismic 
records may be accomplished by exploiting the 
space-time relationship between primary and 
ghost arrivals on a single or multichannel basis. 
An example of the former is provided in Lind- 
sey's (1960) article whereby a simple positive 
feedback loop is employed to effectively cancel 
ghosts from a single near-surface reflector against 
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primaries. In Lindsey's notation" 
apply to each trace is given by 



(I - aH{s)e ") 



filter to 



(1) 



<j = near-surface reflection coefficient. 
//(5) = near-surface filter to account for differ- 
ences in primary and ghost wavelet, 
delay between primary and ghost. 
j= Laplace-transform variable. 

As is evident from equation (I) this method of 
ghost rejection requires specific and detailed infor- 
mation about the ghosting mechanism, which, in 
all but the simplest cases, is difficult to obtain 
from the field records. As Lindsey dennonstrates, 
the autocorrelation function of the input data 
may in principle furnish the parameters a, f/{s), 
and (t) with sufficient accuracy for the feedback 
technique to be employed. This requires, however, 
(1) a high signal-to-noise ratio, (2) wide-band sig- 
nal information, and (3) a simple total seismic 
system wavelet. Any significant amplitude in the 
correlation function due to the primaries them- 
selves at other than zero lag would seriously ham- 
per isolating that part of the correlation due to 
ghosting, and the specification of the parameters 
for equation (1). This will become evident in the 
following. 

The multichannel techniques currently avail- 
able for ghost elimination consist of vertical shot 
stacking (Hammond, 1962), and the use oi dis- 
tributed charges (Musgrave, Ehlert, and Nash, 
1958; Sparks and Silverman, 1953). These, of 
course, differ only in degree and not kind as they 
both represent vertical stacking. The philosophy 
of vertical shot stacking, for ghost elimination, 
finds its basis in the time relationship that exists 
between primaries and ghosts from different shot 
depths. If we have several shots at different 
depths recorded at the same seismometer posi- 
tion, and static correct these traces to a common 
shot depth, the primary reflections will time tie 
and the ghosts, of course, will not. By stacking the 
time-shifted traces, it is clear that the primary 
energy will be enhanced relative to the ght)st 
energy. The stacking may be accomplished in the 
playback center with the individually recorded 
shots, or in the field with delayed charges or dis- 
tributed charges whose detonation velocity is 
matched to the propagation velocity of the sur- 



m 

c 



rounding medium. This is a fl^^^orce ghost- 
cancellation technique which works the ghost re- 
flections against themselves. It is attractive in 
that it does not require knowledge of the ghosting 
mechanism or the primary-ghost relationship. The 
major shortcoming of the conventional shot- 
stacking scheme results from its narrow-band re- 
jection capability when stacking only a few shots. 
This limitation is serious for two reasons; 
(I) There is an industry-wide movement to 
achieve wide-band seismic reflection data for 
fuller utilization of the available signal spectrum; 
therefore, narrow-band techniques are decidedly 
undesirable if alternatives exist; (2) The narrow- 
hand rejection capability makes the shot separa- 
tion a critical parameter; the shot spacing must 
be very nearly { wavelength for the peak fre- 
quency in the spectrum of the ghost wavelet in 
order to dcghost successfully even with narrow- 
band input. 

These considerations lead us to the conclusion 
that the Ideal ghost-elimination processor should 
require: (I) both wide-band ghost rejection and 
signal preservation capabilities, and (2) minimum 
knowledge of the ghosting mechanism and the pri- 
mary-ghost relationship. These requirements can 
l>e well met by allowing for indis idually filtering 
the several shots before stacking as is demon- 
strated in the following. 

METHOD 

This multichannel-processing technique is 
schematically characterized below. 



Where /,(/) is the input trace from the ith shot at 
a given seismometer location, the K/s are filters 
applied to the traces before summation, and g(() 
is the output of the processor. 

The complex transfer functions K,(/) are de- 
signed, on a least-mean-square-error criterion, to 
extract primary events from random noise and 
ghost events, on the basis of the differential lime 
dependence of primaries and ghosts on hole 
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depth. The least-squares-filler design is p d 
in Appendix A, 

Although any number of shots can be stacked 
by this process, two shots usually provide ample 
attenuation of the ghost energy. Therefore, we 
will restrict our attenuation to two-channel prcjc- 
esses. 

In the specification of the signal and noise 
model it is assumed that the signal and noise are 
statistically independent. This assumption insures 
that the resulting filters Ki(/) and VilJ) ^vill not 
require knowledge of the ghosting mechanism tor 
their specification. The fact that the primaries 
and ghosts are in fact correlated does not invali- 
date the process, but rather reflects that wc arc 
not using all the information in the data available 
to us. By utilizing this latter correlation we could 
obtain a somewhat more effective deghostinj; 
processor, Mowever, the difficulty in estimating 
the primary-ghost correlation and the imprac- 
ticality of having to redesign the filters for each 
new ghosting situation presently outweighs any 
potential gain. Before proceeding to examine spe- 
cific examples of the use of the prt)cessor, it is in- 
structive to compare its theoretical response with 
that of the conventional shot-stacking technique. 
The parameter of interest here is the ratio of 
(ghost/primary) amplitude response as a function 
of frequency. Clearly, we desire this to he as small 
as possible consistent with signal preservation and 
random noise rejection at all frequencies. If we 
define the signal and noise for the two shots as 

. = Pit) 

(primary arrivals from shallow shot) 

= no 

(primary arrivals from deep shot) 

and 

.V,(/) = GO) 

(ghost arrivals from shallow shot) 

-V,(0 = GO - At) 

(ghost arrival from deep shot) 

where At is the ^jhost reflection moveout between 
the two shots after static correcting the traces to 
a common depth, and equals twice the difference 
in uphole (2At^) time; then the responses for pri- 



maries and ghosts for the conventional are, 
respectively, 

SM +S,0} = iro) 

.V.(/) + X,{i) = GO) + G{1 - M). (2) 

Transforming equation (2) to the frequency do- 
main we have 

response for primaries 

= 2?(/) 
response for ghosts 

= (7(/)(l -I- (3) 



where 



and 



-f(/) 



-r. 



x{l)e.~'-'fdl 



arc Fourier-transform pairs. 

From equation (3), the transfer function of the 
conventional stack process for primary events is 
the constant 2, and for ghosts ( 1 -f <J~'^'^^') ■ The 
amplitude of the ratio of ghost-to-primary-trans- 
fer functions is clearlv 



I cos rrjAt j , 



(4) 



a well-known result. 

Similarly, the responses for primaries and 
ghosts for the present process, which we call "Op- 
timum Wide-Band Stack" process, are, respec- 
tively, 

response for primaries 

= ?(0(r,(/) + !',(/» 

response for ghosts 

and the amplitude of the ratio of ghost-to-pri- 
mary-transfer functions is 



I K,(/) -t- Y,{f) I 

i YAf) + Y,{f)e-'">^'\ 



(5) 



Expressions (4) and (5) are plotted in Figure 1. 
The ditferencc in the breadth of the reject bands 
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of the two processes is apparent. ^^^^rticular, 

for a AT* of 10 msec, the —20 db reject range for 
the conventional two-shot stack is only 47 lo 5.\ 
cps, whereas the present process provides a 20-80 
cps range. The latter increased the rejection band- 
width by a factor of 10 over the conventional two- 
channel stack. This difference, as well as the gen- 
eral effectiveness of the method, is best illustrated 
in the following deghosting examples using both 
synthetic and field examples. 

RESULTS 

The first example is a highly idealized synthetic 
case, but it serves to put the results of Figure 1 
in the more familiar light of wavelets that might 
be expected to occur on a reflection seismogram. 

Figure 2 shows the idealized cross section in- 
cluding shot and recording geometry used to 
combat ghost reflections on a multichannel basis. 
Two or more shots are recorded at each seismom- 
eter group for slacking purposes. The shot sepa- 
ration is determined by a rough knowledge of the 
velocity in the shot interval, and the center fre- 
quency of the signal spectrum. 



Burg, and Backus 

Figure 3 depicts a possible pr , and ghost 

reflection resulting from the shots A and B re- 
corded at one seismometer group. The traces A 
and B have been time-shifled to a common shot 
depth, thereby lining up primary e\ents. The 
ghost reflections in this case show a 10-msec 
nnoveout or a difference in uphole lime ^Tuh = ^ 
msec. The figure also shows a comparison of the 
conventional stack and the present system tor 
ghost rejection. The former method simply adds 
the two traces with a resulting attenuation of 
ghost amplitude of about —7 db, while the latter 
filters each trace with different filters and then 
stacks, achieving in this instance better than 
— 25 db rejection of the ghost energy. 

The primary and ghost wavelets in this example 
are identical but reversed in polarity. Their am- 
plitude spectrum is given by 




0 < / < 100 



P(f) =0 / > too cps, 




r AT 

Fig. I. Comparison of the present process and conventional two-shot stack response for primary and ghost energy. 
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Fig. 2. Idealized ghost problem in reflection seismology. 
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Fic. 3. Comparison of the present process and conventional stack for ghost elimination utilizing synthetic wavelets. 
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which is peaked at 50 cps. With S^^^ 10 msec 
for the ghost wavelets, it is evident that this spec- 
trum would have its peak frequency centered at 
/pcavAr= .5 in Fij^ure 1. Thus, the results in Fij^- 
ure 3 represent the maKimum rejection possible 
for either process using this wavelet. The advan- 
tages of the wide-hand rejection and signal -preser 
vation characteristics of the process are self-e\ i- 
dent. In addition, no specific knowledge of the 
primary-ghost relationship or the spectrum of tlie 
wavelets, other than the AT = 2ATuh, is required 
to implement the filters Ya and Vit. Fortunately, 
the difference in uphoie time can be obtained ac- 
curately in the field with little ditficulty if reason- 
able care is taken to record an uphoie pulse with 
each shot. Other data which can provide an inde- 
pendent check on the AT" are: (1) accurate knowl- 
edge of the shot separation and vertical velocity 
in that depth increment, (2) correlation functions, 
(.5) difference of arrival times of both primaries 
and ghost from shots at different depths. 

The second example uses field data taken from 
an uphoie survey near Sherman, Texas. 

The spread geometry and details of the uphoie 
survey are shown in Figure 4. Two holes were 
shot from 180 ft to 40 ft in 20-ft increments using 
5-lb charges. The recording array was designed to 
reject the prominent surface waves, based on 
previous noise studies in the area. Only 12 of the 
available 24 channels were used on the surface 
spread; the remaining 12 were employed witli a 
vertical array of seismometers cemented in the 
ground as part of an auxiliary study. Programmed 
gain control and a recording passband of 10-350 
cps were employed. 

The upper half of Figure 5 shows the center 
four traces obtained from SP 2 at the 180-ft and 
160-it shot depths. They have been corrected for 
normal moveout and static corrected to the 
deeper shot. Several strong reflection complexes 
are evident from .2- 4 sec and commencing again 
at about 1.3 sec. The shallow marker at .31 sec, 
and the first part of the reflection complex at 
about 1-4 sec, correlate well between records from 
different shot depths indicating these events are 
primaries (solid circles); however, the remainder 
of the reflections do not have good character cor- 
relation nor do they necessarily time-tie. This is 
particularly evident for the events between .34- 
.40 sec. They lag the primary at .31 sec by less 
than or equal to 27'«ft, or 82 and 76 msec for the 
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180-ft and 160-ft shots, respeclHB .As will be 
demonstrated, these events are ghost reflections 
from the base of the weathering (open circles) and 
surface reflectors. 

The lower half of Figure 5 shows the results of 
conventional and optimum wide-)>and shot stack- 
ing of the 180-ft and 160-ft records. The results 
of the conventional stack indicate fair improve 
ment in ghost reduction over the input traces 
directly above. This is predictable from the two- 
shot-stack response for ghost energy in Figure 1. 
The reflection peak frequency of the records 
varies from about 50 cps for the early events to 
about 30 cps for the later arrivals. The ghost 
moveout for the 2()-ft shot interval is about 6 
msec as can be seen from the records. Therefore, 
the product F^^kAT in Figure I ranges between 
.18 and .30 giving predicted ghost rejection of 
about 3-6 db. The actual achieved rejection ap- 
pears t(» be slightly less than this. The optimum 
wide-band stack results shown in the last 12 
traces indicate that ghost energy has been atten- 
uated of the order of 10-20 db, which is in general 
agreement with the response curves of Figure 1 
for the frequencies and moveout involved. In par- 
ticular, it is noted that the low-frecjuency "leggy" 
portion of the deep reflector at 1.5 sec is, in fact, 
ghost energy. This is unattenuated by the con- 
ventional stack, but is well rejected by the pres- 
ent process. Both shot-stacking processes repre- 
sented in Figure 5 preserve signal (primary) am- 
plitude and waveform, but they differ markedly 
in their rejection capabilities. Furthermore, the 
only information necessary for the implementa- 
tion of the present process was the ghost moveout, 
or twice the difference in uphoie times = 6 msec. 

The very nature of ghosting makes its positive 
identification on the basis of one shot depth al- 
most impossible, even though the ghost reflec- 
tions may be strong. It is felt this is the primary 
reason that the importance of ghost reflections 
has been a somewhat controversial subject in 
reflection seismology. If, however, the seismo- 
gram is viewed as a functi<in of shot depth, the 
primary and ghost energy are clearly separable 
on the basis of their arrival-time relationship. To 
illustrate this, a representative trace from each 
shot depth at SP 2 was gathered onto one record 
to facilitate the shot depth comparison. Rather 
than select a particular trace from each depth, 
however, a velocity filter (Embree, Burg, and 
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Fig. 4. Spread configuration for uphole survey taken near Sherman. Texas. 



Backus, 1963) was applied to each 12- trace record 
to yield a single-trace representation of each shot 
depth. The application of the velocity fitter here 
achieved additional rejection of the horizontally 
travelling shot-generated noise that was not cut 
out by the seismometer groups. Thus, each sinjjle 
trace representation is essentially free of hori- 
zontally traveling energy. What remains is com- 
posed of primaries, multiples, ghosts, and any 
scattered energy with less than two msec per trace 
moveout across the spread. 



The gathered traces are shown in Figure 6 as 
the top eight traces. They have not been static 
corrected to a common depth. Therefore, pri- 
maries will arrive earlier from the deep than from 
the shallow shots while the opposite is true for 
ghosts. The presence of both primary and ghost 
moveout is evident with the formation of the 
characteristic 'T*' pattern due to generically re- 
lated primaries and ghosts. Several of these can 
be seen at record times of .-^4 sec, 1.45 sec, 1.7 
sec, and 1.86 sec. There ran be no doubt on this 
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type of a display that ghost ener^^^^Do^^ pres- 
ent and strong, whereas single depth information 
would not be nearly as revealing. The effect of 
shot depth on character correlation when ghost 
energy is present can best be seen in the reflec- 
tion complex at about 1.45 sec. Notice the dras- 
tic changes in waveform across the 140 feet of shot 
depth resulting from primary-ghost interference. 
Similar horizontal variations also can be expected 
in the primary -ghost relationship for reasons 
previously mentioned. The removal of this ghost 
interference is mandatory before meaningful sig- 
nificance can be attached to reflection waveforms. 
The present process was applied to adjacent pairs 
of uphole traces resulting in the deghosted traces 
immediately below them in Figure 6. The first 
deghosted trace results from "optimum stacking" 
the 180-ft and l60-ft traces, etc. The most strik- 
ing result is that the deghosted traces show only 
primary moveout. The obvious ghost-picks on the 
input traces have completely vanished for visual 
purposes, which implies about —20 db rejection. 
Other significant features are; (1) Several shallow 
reflections have been enhanced. (2) A number of 
weak reflectors in the quiet zone .35-1.2 sec have 
become pickable, (3) Character correlation is ex- 
cellent across the 140 ft of shot depth. The latter 
is particularly noticeable at about 1.4 sec where 
the correlation is equally good for the high-fre- 
quency reflection sequence onset as it is for the 
low-frequency tail. This again illustrates the wide- 
band capability of the processor. The high degree 
of correlation is in no way generated by the proc- 
ess of stacking. In fact, alternate deghosted traces 
have no common data between them. This degree 
of correlation is expected if the following condi- 
tions are fulfilled by the data: 

(1) constant charge size, 

(2) elimination of horizontally traveling en- 
ergy (surface waves, etc.), 

(3) receiver position constant for all shots, 

(4) medium in which the shots occur is homo- 
geneous, 

(5) ghost reflections eliminated. 

Conditions (1) and (3) were accomplished in 
the field, (2) was achieved by the recording array 
and velocity filtering, (4) implies that the same 
wavelet was put into the ground for each shot 
which was substantiated by monitoring the up- 
hole wavelet, and a velocity survey, showing con- 



stant velocity between 30 ft ai^^^D ft, and (5) 

was accomplished by the present process. 

To further evaluate the ghost problem and to 
understand its relationship to the near-surface, 
the process was applied to the data to enhance 
ghost energy and reject primaries. This process- 
ing is shown as the last seven traces in Figure 6. 
The amplitude of the deprimaried traces has been 
boosted by about 6 db to bring them up to the 
level of the input and deghosted traces. Notice 
there is only ghost moveout present. Further- 
more, the ghost reflections appear to be simply 
related to the primaries. For example, the strong 
ghost at about ,35 sec is the mirror image of its 
primary, and has the correct delay for a ghost off 
the base of the weathering. There is a second 
ghost following the aforementioned by 30 msec, 
which would correspond to a ghost ofT the sur- 
face. Thus, it appears that the ghosting mech- 
anism in this area consists of two specular reflec- 
tors, one at the base of the weathering layer and 
the other at the surface. These observations may 
be confirmed in a quantitative manner through 
the use of both correlation functions, and the esti- 
mation of the impulse response of the near-surface 
obtained from computing the optimum linear 
operator for transforming the deghosted traces 
back into the original ghosted traces. 

GHOST ANALYSIS USING CORRELATION FUNCTIONS 

The use of the autocorrelation function, and its 
interpretation for a single ghosting horizon has 
been discussed by Lindsey (1960). His analysis 
can be readily extended to more complex ghosting 
situations. Of particular interest here is the pre- 
dicted autocorrelation function for a double ghost 
mechanism. 

Let the trace represented by primary energy 
only be p{t). The ghosted trace resulting from 
two near-surface reflectors (above the shot) is ap- 
proximately given by 

R(0 = p{t) - ctip(( - n) - ppit - rj) (6) 

T| — two-way traveltime from shot to the first 

ghost generator. 
Ti = two-way traveltime from shot to the sec- 
ond ghost generator, 
ai — reflection coefficient of the first reflector, 
and 

^ = apparent reflection coefficient of the sec- 
ond reflector. 
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Fio. 7. Impulse autocorrelation function for primary plus double-ghost problem. 



The true reflection coefficient of the second re- 
flector, art, is related to ^ by at = /C3/( I — Ofi*). In 
the above, multiple reflections between the two 
reflectors have been neglected. 

The correlation function 0(r) for the ghosted 
trace is given by 

I r ^* 

<t>p^Ar) = — I g{t)g{t + T)dt. 

A — TvJ r, 

Substituting in the above for g{t) given by equa- 
tion (6), expanding and collecting like terms, the 
ghosted autocorrelation function may he written 
as a sum of seven terms: 

^p^.ir) = (1 -h a,» -h ^')<^p(r) 

-f ai^0p(r — r2 + ti) 

+ ct\^4>p{r + Tj — Ti) 

— cc\<t>p{T — ri) — ai<^p(r + ti) 

- /S<^p(r - T,) - 0<t>p{r -f r^) (7) 

where 

<^p(0 = — piOpO + r)dt 

is the autocorrelation function of the primary re- 
flection sequence. Expression (7) may also be 
written as the convolution of 0p{r) with a se- 
quence of impulses as 



{t)h{r - l)dl. 



with 

A(t) = (I + ar + J3')5(r) + a^S{T - r, + r.) 
+ ai/J5(r + Tt — Ti) — QfiS(r — ri) 

- a,5(r + T,) - /34(r - rj) 

- ;3«(r - n). (8) 

It is more instructive to examine A(r), the impulse 
autocorrelation, than <^p^.o(r) since the latter con- 
tains the smearing effects of the reflection wavelet. 

The impulse autocorrelation function h{r) for 
primary and double-ghost situation is shown in 
Figure 7. Note that if j3-*0, the autocorrelation 
function reduces to that given by Lindsey (1960) 
after convolution with the primary autocorrela- 
tion function 0p(r). The addition of a second 
ghost reflector results in two additional features 
in the structure of the correlation function as 
seen in Figure 7. The negative spike of amplitude 
/3 at ±rt sec results from the correlation of pri- 
mary and ghost reflections separated in time by 
T2 sec, that is, the second set of ghosts; while the 
positive feature of amplitude ai/3 at ±(ri— ti) 
seconds is due to the correlation between the two 
ghosts themselves. The degree to which these 
features may be seen on an autocorrelation func- 
tion depends, of course, upon the resolving power 
of the primary correlation </>p(r). It the latter is 
very nearly a delta function itself, then our pic- 
ture in Figure 7 would remain essentially un- 
changed. If, however, 0p(r) is slowly decaying 
away from zero lag because of reverberations or 
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narrow-banti signal information, tWBKir impulse 
representation would be severely distorted after 
convolution, and our ability to identify particular 
wiggles in the correlation function with ghosting 
would be practically nil. 

The autocorrelation functions were computed 
for the input traces in Figure 6 and are shown in 
Figure 8 arranged in order of decreasing hole 
depth. Only the positive half of the correlation 



function is shown, it being sym^^PTc about zero 
lag or T — 0. Twice the uphole time, 2T,j^, as 
determined from the uphole pulse is indicated by 
an arrow below each correlation function. There 
is generally good agreement between this time 
and a trough in the correlation functions pre- 
dicted for the primary-surface ghost correlation. 
A second stronger trough can be seen preceding 
the surface ghost trough by about 30 msec, and 
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Fic. 8. AutocorrcUtion functions of uphole traces. 



BNSDCXitD: <XP 2 06781 9A I > 



Technique for Elltn ^natlon of Ghost Arrivals 

1 1 r 




795 



-f ( LAGS) mmmc 



Fig. 9. Autocorrelation functions of deghostcd traces. 



paralleling the latter's moveout between shot 
depths. This trough is identified with a base-of- 
the-weathering ghost reflection, which is the most 
prominent ghost on the records. There is also a 
f>ositive correlation at about 30 msec independent 
of depth, that is, predicted by Figure 7 for two 
ghosts separated by 30 msec. 

The correlation functions of the deghosted 
traces of Figure 6 are shown in Figure 9. The two 



troughs associated with ghosting which move 
across the correlations in Figure 8 are no longer 
evident- In fact, the residual wiggles show good 
reproducibility from shot depth to shot depth 
with the exception of the shallower shots. They 
represent the structure in the autocorrelation of 
the effective seismic wavelet, <^>p(t), which is inde- 
pendent of shot depth for reasonable depths. 
Several additional features become clear with 
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the before and after correlations ^|^fures 8 and 

9, The l2()-ft ghosted correlation function shows 
anomalous lack of correlation for r greater than 
30 msec. At this depth, however, the ghost 
troughs coincide with peaks in the primary corre- 
lation ^p{r) causing destructive interference. The 
primary autocorrelation functions in Figure 9 do 
have an intrinsic peak at about iO-msec lag, how- 
ever, this peak is reduced from Figure 8 indicat- 
ing that the positive correlation of the two ghosts 
contributed to its amplitude. 

It becomes evident what difficulties might be 
encountered in attempting to use amplitude infor- 
mation from the ghosted correlation functions to 
estimate reflection coefficients and near-surface 
filtering effects on the ghost reflections. Not hav- 
ing to specify these parameters makes shot-stack- 
ing a particularly attractive deghosting tech- 
nique. 

ESTIMATE OF NEAR-SURFACE REFLECTION RESPONSE 

Once we have obtained a ghost-free estimate of 
the subsurface reflection sequence, we can work 
backwards and achieve a quantitative descrip- 
tion of the near-surface reflection complex. This is 
accomplished by computing the operator or filter 
which transforms a deghosted trace into one of 
the ghosted traces from whence it came. The im- 
pulse response or time-domain representation of 
this operator should embody the effects of the 
near-surface on the reflection seismograms, e.g., 
the cause of the ghosting. 

Mathematically, we seek the operator /i(r) 
which will transform a ghost-free trace p{t) back 
into a ghosted trace git) or which satisfies the 
relation, 



«(0 



P(r)hU - r)dr. 



(9) 



At this point the reader may wonder if he is not 
caught in a vicious circle of first removing ghosts 
and then replacing them in the data. To expel 
any such fears a word of explanation is in order. 
Recall that the present process by design only 
requires knowledge of the difference in uphole 
times for its implementation. No other informa- 
tion about the ghosts is required such as: reflec- 
tion coefiFicients, primary -ghost time lag, near- 
surface filtering effects, or details of the ghosting 
mechanism. Therefore, these parameters which 
are interesting from the standpoint of studying 



the near-surface and its effec^^^he reflection 
seismogram do not automatically drop out of the 
processing. They may be inferred, of course, by 
visually examining the traces before and after 
deghosting, and computing correlation functions 
and the like. The integral e(iuation (9) however, 
provides a quantitative measure of these near- 
surface parameters in terms of the unknown func- 
tion h(t). 

The least-squares estimate of ft{i) may be ob- 
tained from g(0 and the deghosted estimate of the 
primary reflection sequence pU). That is, we seek 
the filter h(t) which minimizes the mean-square 
error, 

= ^ f kit) - Ho * pO)Vdi (10) 



where * denotes convolution. Solutions of equa- 
tion (1(1) for discrete data are given by Lcvinson, 
1947, and more recently by Levin, 1960, and 
Robinson, 196.i. 

In matrix form, the least-squares estimate of 
Ht) is, 



♦p(0) M2) 

*pO) MO) 

4>pIn) 



A(2) 



(0) j,A(iV)j 
«r,(0) 1 
*P»(1) 
«P,(2) 



(11) 



«P.(/V), 



where <^p(i) is the tth lag autocorrelation oip(i)^ 
and <t>po(i) is the tth lag crosscorrelation of ^(0, 
and g{t). The special symmetry properties of the 
correlation matrix in (11) allow an iterative solu- 
tion for the filter w^eights A(t). That is, the solu- 
tion vector can be obtained by operations 
on the previous vector (A(t — 1)]^ without actually 
inverting the »X« correlation matrix, as demon- 
strated by Levinson, 1947. 

Figure 10 shows the operator calculated 
for several different shot depths at SP 2. 

At each depth represented, the operator was 
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computed from a deghosted trace and the gh 
trace corresponding to that depth. The spike in 
the operators at / = 0 simply reproduces the pri- 
mary events when convolved with ^(0- The first 
strong negative trough produces the base-of- 
weathering ghosts, and the second shallow, 
broad trough correlating with twice the uphole 
time gives rise to the surface ghost. As also ob- 
served on the autocorrelation functions, this 
double set of troughs is separated by 30 msec. 
The base-of-the-weathering trough remains con- 
stant up the hole in both amplitude and charac- 
ter, whereas the shallow surface ghost trough 
exhibits sizeable waveform fluctuations from 
shot to shot. This may he indicative of the noise 
level of the data. 

The reflection coefficients cti and 0 defined in 
equation (6) were estimated from the operators 
in Figure 10 by taking the ratio of ghost trough to 
primary peak at / = 0. .Admittedly, this may be a 
crude approximation to the true reflection coeffi- 
cients in view of the difference in frequency con- 
tent between the base-of-weathering and surface 
ghost trough. Clearly and oti are frequency de- 
pendent, or at least our simple mode! of two 
sharp ghosting horizons forces them to be. Since 
we have not allowed for frequency dependent 
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Fig. 10. Ncar surfacc response estimate ob- 
tained for several shot depths. 




Fic;. 11. Near-surface response estimate obtained 
from transforming the l80-ft deghosted trace into the 
1 20- ft original trace. 



absorption within the weathering layer, any 
effects due to absorption are necessarily lumped 
in the parameter 0. Keeping these facts in mind, 
the average reflection coefficients obtained from 
the near-surface responses in Figure 10 arc: 

ai = .47, 

^ = .28, 

giving for the surface reflection coefficient 



1 - ai' 



.28 
.78 



= .36, 



These estimates provide a lower limit for the 
reflection coefficients because of the imperfect 
rejection of ghost energy and because of noise cor- 
relation between deghosted and original traces. 

To reduce the effect of residual ghost and 
noise correlation, the optimum filter to transform 
the 180-ft deghosted trace into the 120-ft original 
trace was computed. The filter, limited to about 
70 cps on the high side, is shown in Figure 11. 
The reduction of "primary'* correlation relative 
to ghost correlation is apparent. Crudely com- 
pensating for the low-frequency drift in the oper- 
ator, the values. 



.4 
.6 



are obtained. This yields a model indicated in 
Figure 12. 
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Fic. 12. Two-interface mottel for near-surface. 




Fig. 13. Crosscorrelation functions computed between different shotpoints for the 
original vclocily-fiUcred uphole traces. 
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The fact that the effective reflecti<)n cot^^Hlt 
of the surface is less than one (i.e. ^, = .6) is pre- 
sumably due to energy dissipation in the weather- 
ing layer. This is further supported by the appar- 
ent lack of high frequency (above 50 cps) in the 
weathering ghost. 

The simple two-interface discrete model shown 
in Figure 12 yields a result that the total energy 
in the ghost complex should be about 54 percent 
of the total initially down-traveling wave, 
(.6^-f .V+.l5'-|-.06« = .54). The adequacy of this 
simple model can be further tested by examina- 
tion of the crosscorrelation function computed 
between different shot depths for the original 
velocity filtered traces, shown in Figure 13. As 
exemplified in the crosscorrelation between the 
180-ft shot and the 40 ft shot, all of the primary 
energy should correlate at a lag r— —ATj^, and 
all of the initially up-traveling energy should cor- 
relate at a lag r = -\~AT^h, no matter how complex 
the ghosting mechanism is. The 180-160 ft cor- 
relation is uninterpretable because of the small 
ATuA involved. The 20-ft shot is in the weathering 
and yiehls a very unreliable correlation. The two 
center correlations in Figure \^ suggest that the 
initially up-going energy contributes about 80 
percent as much energy to the reflection record 
as the initially down-going energy. This com- 
pares with a .S4 percent figure accounted for by 
anal^'sis on the basis of the simple two-reflector 
model in Figure 12. The 80 percent figure sug- 
gests that about 20 percent oi the initially up- 
going energy in the reflection-record frequency 
range is dissipated in the weathering layer or scat- 
tered by undulation in the boundaries of the up- 
per reflectors. 

The major portion of the analysis was focused 
on SP 2. However, the SP 1 data yielded similar 
findings, with the exception that the ghosting 
appeared somewhat more complex. This was evi- 
dent on the uphole traces as well as the autocor- 
relation functions, which were less interpretable 
from a ghosting standpoint than those in Figure 
8 from SP 2. After deghosting, the SP 1 traces 
correlated very well with the SP 2 deghosted 
traces. 

The estimate of the near-surface reflection re- 
sponse at SP I obtained from transforming de- 
ghosted into ghosted traces revealed differences in 
the near-surface response between SP 1 and SP 2. 
Figure 14 shows a comparison of the near surface 
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Fig. 14. Comparison of near-surface res[)on5e estimates 
obtained at SP 1 and SP 2 for the 180-ft shot. 

response estimates obtained at the two shot- 
points lor the 180-ft shot depth. The portions of 
the operators relating to primaries are almost 
identical. Twice the uph<)ie times indicated by .S" 
are significantly tlifferent at the two shotpoints, 
yet the base-of-the-\veathcring reflector remains 
essentially fi.xed, therefore, the traveltirne through 
the weathering at SP 1 is about 7 msec shorter 
than at SP 2. The base-of-the-weathering ghost 
trough at SP I is of comparable strength to the 
SP 2 trough, but docs not possess the symmetry 
of the latter. Furthermore, 2T,j^ = S at SP I does 
not correspond to a well-defined ghost trough ; and 
finally an additional strong positive feature ap- 
pears at the first reverbt.-ration time within the 
weathering layer. Therefore, it is likrly that the 
total interval from about 50 to 90 msec at SP 1 
represents ghosting, the details ol which are not 
resolved as simple troughs or peaks but rather 
form a complex interference. 

CONCLUSIONS 

A new deghosting technique is presented which 
utilizes two shots at different depths. It differs 
from conventional stacking in that separate filters 
are applied to the two shots before stacking. The 
filters are designed in an optimum manner to en- 
hance initially down-traveling energy and reject 
initially up-traveling energy. The resulting proc- 
ess has wide-band ghost reject and signal-preser- 
vation characteristics. \o specific knowledge of 
the ghost mechanism is required for its imple- 
mentation other than, (1) determining the differ- 
ence in upholc times, (2) assuming a spectral 
content fi>r signal and noise, usually ''white'* over 
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Fig. is. Left half of the figure shows the sampled impulse response of the deghosting fillers, and the right half shows 
the impulse response for primaries and ghosts obtained by suitably adding the filter responses. 



the frequency range of interest, and (3) insuring 
that both shots are below the ghosting mech- 
anism. The method will handle ghosting from a 
reflection complex or ghosts from a single reflec- 
tor equally well. 

The degree of primary and ghost energy sepa- 
ration achievable with the present process permits 
the quantitative investigation of the near-surface 
reflection response by computing the operator 
which transforms the deghosted into the ghosted 
trace. 
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APPENDIX A 
LEAST-SQUARES TWO-CHANNEL PILTER DESIGN 

The classical least-squares technique is applied 
to two-channel sampled time functions for opti- 
mum-filter design by an obvious extension of 
Levinson's (1947) single-channel treatment. 

Consider the two-channel time functions fi(t) 
and /t(0 to be represented by their sampled 
values at points l = kh, where h is the sample pe- 
riod and k is an index. Let the sampled functions 
j\{kh) and /zikh) be regarded as the sequences bit 
and Ajfc representing signal and noise, that is: 

l>2k — Stk ~h A^2*> 
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where 



Raik) = lim 



Sik = signal on channel I, 
= signal on channel 2, 
A'u = noise on channel 1, 
Ntk — noise on channel 2, 

We wish to determine the pair of linear filters 
which, when convolved, respectively, with the in- 
put sequences but and bu and summed, will give ^[q^ (A-2) we have: 
an output as similar as pK>ssibIe to the desired sig- 
nal on one of the input channels^ say 5ut. The fiJ- 
ter sequences Ku and Yu are designed such that 
the total power in the error signal 



AT-. 2.V -h 1 .11^ 
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and the crosscorrelalion function between two 
sampled time functions as 

R^{k) = lim \ E (a.b.^) 
jv-.- 2 A' + 1 ,__.v 

= /?,»(- k). 

Substituting the correlation functions into equa- 



I = /?,,(()) - 2 Z Ki„/?6..,(n + p) 

M 

- 2 ytnR6,.,{n -f p) 



(A-l) 



is minimized. The power in the error signal is 
given by the sum of the squares of the £fc's, or 
slated in formula 



n*0 m — 0 
Si Af 



(A-3) 



lim 



AT^- IN + 1 kZ-!, 

1 



= lim 



jv^. 2yv 4- 1 *r_Y 



Z (5u^p)^ - 2 £ Fu lim — ^ _ 

y-- 2/V + 1 jt — N 



2 Z Hm 



1 



2/V + I iir^v 



+ 2 Z Z Yrr.Y^^ lim — — - 
„-o«.o 2N + 1 

Af Af 



E b\k-nb2 
k .V 



.V 

i« lim E Au-n^IA- 

IN -h 1 .n^v 



W Af 



2„ !im — — 7 Z *2*-n*2Jt-m, 

iv-*- 2A^ + 1 .V 



(A-2) 



which is to he minimized with respect to the Ku's 
and Kat^s. 

Equation (A-2) may be simplified by introduc- 
ing the correlation functions. 

The autocorrelation function for a sampled 
time function is defined as, 



For / to he a minimum with respect to the Ku's 
and Kui*s we must have 

d[ 

A = 0, I, 2 • • . Af 



SYik 
dl 



= 0 
= 0 
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dvZ 



- 2Rt,.X^ + p) 

M 

+ 2 E I^2«/26,6,(m - k) 
+ 2 E Vi„R,,{k - «) = 0 

Wl-O 
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This system of equations may be written in 
matrix form as 



ro Tx ft 



IK,/ J 



7 m 



(A-5) 



where each of the r/s, K^'s, and Tj's are sub- 
matrices defined as follows 



Notice also that r^-k-rk-^. 

An iterative solution to (A-5) is also possible 
(Robinson, 196.5) as in the single-channel case. 
The parameter p in Hi and equation (A-1) relates 
.to whether or not one ficsires to do filtering, filter- 
ing and prediction, or interpolation. If /> = 0, we 
attempt to extract .SV from knowledge of the 6's up 
to time t Ilk. Ti p >0, wt attempt to predict 5*.^p 
ahead in time hy l^hp sec, and if <0, we use 
values of the i's in future time to interpolate the 
value of Sk-p- 



Up until now the derivatio^^^^ been com- 
pletely general without regard to "the deghosting 
problem. In order to design the deghosting filters 
on a theoretical model the various correlation 
functions in (A-4) must be specified. 

The signal and noise description for the ghost- 
reflection problem on a two-channel basis may 
be adequately specified as: 

Shallow shot 



hit) 

Deep shot 



where 

P\.2{t) — pT[m3,xy reflection sequence on chan- 
nel 1, 2, 

^7j j(f) = ghost reflection sequence on channel 
1, 2, 

N\.2{t) = random noise on channel 1, 2. 

If the traces /l(0 and /^{t) are static-shifted 
by the difference in uphole time=drtjn then 
Py{t) = P2(t)^P{t). and G.(0 = G(/), G^{i) 
= G{t-2^Tuh), based on the assumption that 
both shots are in the same medium below the 
ghosting horizon. We further assume that the ran- 
dom noise is white with power density level No, 
and is uncorrelated with primaries or ghosts. The 
final assumption concerns the crosscorrelation of 
primaries and ghosts. Of course, physically they 
are correlated as Figure 8 shows, however, the 
measurement of this correlation is difficult even 
with good data. Furthermore, the primary ghost 
correlation can be expected to change drastically 
in short horizontal distances due to variable near- 
surface conditions as well as variations in shot 
depth. From the practical consideration of filter 
design we neglect the primary ghost correlation 
and accept the slight loss in filter effectiveness. 

With these assumptions, the various auto- and 
crossrorrelations become in the continuous case: 





= Mr) + ■+- -Vo^Cr), 




= *p(0 + <t>,(.r) 4- A'u«(r), 


Rini'ti'^) 


= Mr) + - 2A2'„^), 


RbtbX^) 


= M-r) + <pAr + 2AT^k), 




= Mr), 




= Mr)- 



(A-6) 
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In order that the filters have a wide-band 
bility for separatins primaries an<J ghosts \vc as- 
sume that the spectral content of the primaries 
and ghosts is white over the frequency range of 
interest. That is, we are not introducing any fre- 
quency difference between signal and noise, hut 
rather exploiting only the spatial correlation dif- 
ferences. Thus, the correlations become 
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0,(t) = G5(r), 



(A.7) 



where P and G are the signal and ghost-power 
levels and 5(r) is the delta function. 

Figure 15 shows the impulse response of the 
filters Kim and Kj^ designed by means of equa- 
tion (A-5), using correlations defined by (A 6) 
and (A-7) for the following parameters: 

f* = 1.0 primary power level, 

G — i.O ghost power level, 

No == -1 random noise power level, 

2ATuh = 1.0 times the sample period = h, 

M ~ 20 number of filter points minus one, 

/» = - 10. 

The figure also shows the impulse response for 
primaries obtained by adding Kim4-K>„, and the 
impulse response for ghosts obtained by shifting 
and adding Vim-\- Y-im-i- Note that the response 
for primaries is very nearly a delta function with 
maximum side-lobe levels of less than 10 percent 
of the main peak, whereas the impulse response 
for ghosts is uniformly small. 

The time-domain formulation of the least- 
stjuare filter problem presented here is clearly 
superior to a frequency-domain formulation ii 
the resulting filters are to be used as finite digital 
convolution operators. The reason being that the 
formrr provides an optimum finite length oper- 
ator, whereas the latter yields an optimum ampli- 
tude an<l phase spectrum which must be Fourier 
transformed and in general truncated in the time 
domain. In the limit of infinite-length operators, 
the time and frequency (hjmain results are identi- 
cal. Inasmuch as filter or operator length is an 
important economic factor in the implementation 



of this and other related techniques it isl^^ble 
to achieve the maximum effectiveness of the filter 
with the least number of points. The time-domain 
solution insures that each of the A^-filler points 
used is optimum. In addition, by examining the 
mean-square error as a function of number of filter 
points, it is possible in most problems to select a 
minimum filter length beyond which the mcan- 
st^uare error essentially levels off. That is, the 
addition ol more points to the filter does not sig- 
nificantly improve its performance in the mean- 
square-error sense. 

The frequency domain on the other hand, if not 
best suited for digital filter design, does provide 
physical insight and interpretation that is lack- 
ing in equations (A-5) and (A-6). 

The frequency-domain solution to multiple 
time series least-mean-square-error filter prob- 
lems is contained in Wiener's (1^49) original 
work. Certain important simplifications in the 
mathematics result when both the past and future 
of the time series are available for filter design, 
as is the case of magnetic-tape-recorded seismic 
data. The solution for the optimum filters Vjij) 
operating on both past and future is (Spieker, 
Burg, Backus, and Strickland, 196t; Foster and 
Sengbush, 1964; and Burg. 1964) given by the 
matrix e(juation, 

[Soif) + .V,,(0](K,(/)] = [5a.(/)] 
= 1. 2, • - • , -V, 
> = 1, 2, . . • , .V, (A-8) 

where S^,(f) and N,j(f) are the cross spectral den- 
sities between channels i and j for signal and 
noise, respectively. The column vector [.^*/(7)] 
consists of cross spectral densities between signal 
in channel t and the desired signal to be estimated 
at channel k. 

The Vjif) are, of course, the unknown filters to 
be applied to eacli channel in the multichannel 
signal-extraction process. Fcjuation (A-8) implies 
that signal and noise are uncorrelatcd. An addi- 
tional trrm in the equation would allow for the 
latter if it exists. 

Using the si^^nal and noise model for the ghost- 
ing problem defined h\' i:t[uation (.^ 6), and trans- 
forming to the frecjuoncy domain, we have a 2X2 
matrix inversion problem lor the filters Viif) and 
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*,(/) + *.(/) + No^mM) + 



or, Lapn«r, Burg, and Backut 



*p(/) + 



+ *,(/) + No J Lf,(/)J U,(/)J 



(A-9) 



where 



The solutions of (A-9) are, 



ration. The denominator of (A-II) which is com- 
mon to both Ki(/) and Kt(/) is essentially a com- 
pensation filter to undo the distortion introduced 
in the primaries as a consequence of ghost attenu- 
ation. 

The response for ghosts may be obtained iti a 
similar fashion by adding and 
giving, 

[See A- 13 ai top of next page] 



7« + rc(l - 



7r(7« + 2yc + 2) + 27.(1 - cos IwflAT^f) 
Y^if) = Yi*{J) * = complex conjugate. 



(A~10) 



The 7*s 
that is, 

7r 



are the noise- to-signal-power ratios, 



7c = 



random-noise- to-signal-power ratio 



= coherent-noise-to-signal-power ratio. 

As previously mentioned, the response for pri- 
maries is obtained by adding and Yt{J) 
which yields, 



The latter is small at all frequencies except 
the set /n, at which the primary and ghost re- 
sp<jnses are equal as evidenced by (A- 1 3). At 
/„^m/4AruA, m=== 1, 3, 5, 7, ■ - equation (A-13) 
has a zero irrespective of the value of yn and 7c. 
At these frequencies the ghost energy is 180" out 
of phase between the two shots and cancellation 
is effected by simply adding the two traces. This 
frequency corresponds to the node in the response 
curves presented in Figure 1. The latter were ob- 
tained from equations (A ll) and (A-13) with 
7ft = .1 and y^ — 1.0, 

It would seem that by allowing the random- 
noise-to-signal-power ratio 7/f to become arbi- 



2yH 4- 2-Y.(l - cos lirflhT^^) 
•y«(T« 2y, + 2) 4- 2y.(1 - cos 2T/2Ar.O 



(A-11) 



For 7ie«7« the response for primary energy is 
very nearly equal to unity for all frequencies ex- 
cept at and near the set /n = n/2ATuJt, n = 0, 
1, 2 - • • where the response becomes, 



K.(/„) -I- = 



(>« -1- 27. + 2) 



(A-12) 



The /„ are alias frequencies at which the stack 
(either conventional or the present process) can- 
not distinguish between primary and ghost 
energy. The associated wavelengths \n—2d/n 
are submulliples of twice the vertical shot sepa- 



trarily small one could effect perfect undislorted 
separation of primaries and ghosts- That is, equa- 
tions (A-11) and (A-13) tend to unity and zero, re- 
spectively, as 7«— »0. The difficulty, of course, is 
that the amplitude density of the filters (A- 10) be- 
comes infinite at /^/^ for 7« = 0. This implies that 
the impulse response of the filters is not square 
integrable over the time interval (— «). In- 
deed, the impulse response of the filters for7«^0 
can be obtained by expanding the periodic fre- 
quency spectra (A-IO) in a Fourier series and 
equating the coefficients Om to delta function 
amplitudes at t=-m2ATuh. That is, 
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2yR cos Itt/ATuh 



yniyR + 27c + 2) + 2tc(1 - cos 2x/2Ar,0 



.t: 



805 

-13) 



A,(0 = Z ''im^C/ - m2Ar„0 (A-U) 



where 



+ i form < 0 
f or m > 0 
-i form < 0 
+ i form > 0 



clearly, 



= response for primaries = 



= response for ghosts = 0. 

These filters do not decay with increasing or 
decreasing tin^e, and are obviously not practical 
for digital application. The role of random noise 
in (A- 10) clearly provides stability and conver- 
gence for the impulse responses such that they 
may be truncated after a finite number of samples 
without degradation of the fillers. The time- 
domain formulation circumvents the problem of 
having to truncate the filters at all by solving 
(A-5) for optimum A^-point filters. The latter 
will, of course, depend upon 7^ as well as the 
other parameters used, however, they are the best 
A^-point operators in the mean-square-error 
sense. 
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